o 
o 

1—5 

m 






X 



Stochastic and Discrete Time Models of Long-Range Turbulent 
Transport in the Scrape-Off Layer 



D. Volchenkov * R. Lima 



Centre de Physique Theorique, CNRS, Luminy Case 907, 



^ ■ 13288, Marseille CEDEX 09 France 

a , 
'"q^! E-Mail: dima427@yahoo.com, lima@cpt.univ-mrs.fr 

O ; February 2, 2008 

C/2 
>>. 

^, 

'~^' Abstract 

Two dimensional stochastic time model of scrape-ofF layer (SOL) turbulent transport is studied. Insta- 
bility arisen in the system with respect to the stochastic perturbations of both either density or vorticity 
reveals itself in the strong outward bursts of particle density propagating ballistically across the SOL. The 
C^ , stability and possible stabilization of the cross-field turbulent system depend very much upon the reciprocal 

/^ I correlation time between density and vorticity fluctuations. 

^H ' Pdf of the particle flux for the large magnitudes of flux events can be modelled with a simple discrete 

^^ ■ time toy model of random walks concluding at a boundary. The spectra of wandering times feature the pdf 

Y\ ■ of particle flux in the model and qualitatively reproduce the experimental statistics of transport events. 
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1 Introduction 

Turbulence stabilization in plasma close to the wall blanket of the ITER divertor is the 
important technical problem determining the performance of the next step device. Long range 
transport in the scrape-off layer (SOL) provokes the plasma-wall interactions in areas that are 
not designed for this purpose. Evidence of the strong outward bursts of particle density prop- 
agating ballistically with rather high velocities far beyond the e-folding length in the SOL has 
been observed recently in several experiments (see PQ-jSI) and in the numerical simulations ^. 
These events do not appear to fit into the standard view of diffusive transport: the probability 
distribution function (pdf) of the particle flux departs from the Gaussian distribution forming 
a long tail which dominates at high positive flux of particles |4| . 

Theoretical investigations of the reported phenomena remain an important task. In the 
present paper, we consider a variety of two dimensional fluid models based on the interchange 
instability in plasma studied in jSj-jH] and discussed recently in |4j exerted to the Gaussian 
distributed external random forces to get an insight into the properties of turbulent transport 
in the cross-field system. 

Neglecting for the dissipation processes in plasma under the constant temperatures Tf, ^ Ti, 
this problem is reduced to the interactions between the normalized particle density field n{x, y, t) 
and the normalized vorticity field w{x,y,t) related to the electric potential field (l){x,y,t), 

dtn = [n, 0], dtw = [w, 0] - gdylogn, w = A^0, g>0, (1) 

defined in the 2D plane transversal to e^, the direction of axial magnetic field Bq. In (^Q), a; 
and y are the normalized radial and poloidal coordinates respectively. The Poisson's brackets 
are defined by [p, q] = dxpdyq — dxqdyp. 

When g = 0, the equations (jT)) describe the 2D-rotations of the density and vorticity gradi- 
ents around the cross-field drift v = —c/Bq V0 x e^, in which V = (9^, dy). Their laminar solu- 
tions (with ty = 0) are given by any spatially homogeneous electric potential (f) = V'i(^) aiid any 
stationary particle density distribution n = ip2{x, y). Other configurations satisfying ((H) at (7 = 
are characterized by the radially symmetric stationary vorticity fields w = d'^(f){x) with the elec- 
tric potentials invariant with respect to the Galilean transformation 0(x) -^ (f){x)+x(fi{t)+if2{t) 
where the parameters of transformations v^i,2(^) are the arbitrary integrable functions of time 
decaying at t ^ —00. The relevant density configurations n = (^^ (x,y — J_^Vy(x,t') dt') 
have the form of profile-preserving waves convected in the poloidal direction by the poloidal 
cross field drift Vy{x, t). The poloidal component of cross field drift itself remains invariant with 
respect to the Galilean transformation v^^ ^ v^^ -|- f2{t), while its radial component v^ = 0. 



Configurations that satisfy (0) for g > have the Boltzmann density distribution of particles 

in the poloidal direction. In particular, those solutions compatible with the Galilean symmetry 

discussed above (with Vj^ ^ 0) are the solitons (solitary waves) of density convected by the 

poloidal electric drift, 

1 

ncxexp — 

9T{x,y) 

where T{x, y) is an arbitrary function twice integrable over its domain. In addition to them, for 

g > 0, Eq. (^ allows for the radially homogeneous configurations dxn = 0, w = f/(y)mod27r 

with Vy = which do not fit into the Galilean symmetry, these are the steady waves. 



y- Vy{x,y,t')dt' 
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n oc exp — / U(y') dy' . 

9 Jo 



The latter solution does not possess a reference angle and can be considered as an infinitely 
degenerated state of the system since the relevant configurations {n, w} can be made equal 
at any number of points by the appropriate choice of U: U{yi) = U{y2) . . . = U{yn), and 
fo^ U{y') dy' = J^^ U{y') dy' . . . = f^" U{y') dy' . For instance, it can be represented by the 
periodic lattice potential controlled by the spokes of high particle density radiating from the 
center. With two concurrent symmetries there can occur either the frustration of one of them 
or the vanishing of both with the consequent appearance of a complicated dynamic picture 
that is most likely stochastic. The latter case corresponds to a maximally symmetric motion 
resulting from the destruction of unperturbed symmetries [Zj . In particular, instability in the 
system (P) occurs either with respect to any small perturbation of density or vorticity. 

Accounting for the dissipation processes in plasma smears the picture, so that the small scale 
fluctuations would acquire stability. In the present paper, we demonstrate that the small scales 
fluctuations can be stable provided there exist the reciprocal correlations between the stochastic 
sources of density and vorticity in the dynamical equations. The large scale stability of a 
fluctuation can be characterized by the order parameter ^ = \ky\ j [k^ + fc^) in the momentum 
space where k^ and ky are the radial and poloidal components of momenta respectively. For 
the uncorrelated random forces (under the white noise assumption), a fluctuation with ^ > 
is unstable with respect to the large scale asymptote in the stochastic problem. 

The accounting for the convection of particles by the random vortexes introduces a finite 
reciprocal correlation time Tc( |r — r'| ) between the density and vorticity random forces. Then 
there exists the critical scale ^c, in the stochastic model, such that a fluctuation with ^ < ^c 
vanishes with time, but its amplitude grows up unboundedly with time, for ^ > ^c- 

Biasing of wall components can locally modify turbulent transport and is considered to be 
beneficial if one aims to insulate the tokamak main chambers from the bursts of density jlj. 



Indeed, the generation of a uniform electric drift in the poloidal direction, v^ — > v^ — V, would 
frustrate one of the symmetries in (^ reestablishing the Galilean invariance in the system. 
For instance, those configurations characterized by the trivial poloidal component of electric 
drift Vj^ = would be eradicated. In the present paper, we have investigated the problem 
of turbulence stabilization close to the divertor wall in the first order of perturbation theory 
and shown that there exists a critical value \Vc\ < oo of the poloidal electric drift which 
would suppress the large scale instability in the stochastic system with the correlated statistics 
of random forces, ^c > 0. However, for the uncorrelated random sources in the stochastic 
problem, ^c = and \Vc\ ^ oo as k -^ 0. 

Correlations between the unstable fluctuations of density at different points are described 
by the advanced Green's functions which are trivial for t > 0. In particular, these func- 
tions determine the concentration profile of the unstable fluctuations of density which increases 
steeply toward the wall. The size of such fluctuations grows linearly with time. In this case, 
the statistics of the transport events responsible for the long tail of the flux pdf is featured 
by the distribution of the characteristic wandering times of growing blobs convected by the 
highly irregular turbulent flow in the close proximity of the divertor wall. In our model, we 
have replaced this complicated dynamics with the one dimensional (the radial symmetry is 
implied in the problem) discrete time random walks. Such a discrete time model would have 
another interpretation: the advanced Green's function is a kernel of an integral equation which 
relates the amplitudes of the growing fluctuations apart from the wall with those on the wall, 
in the stochastic dynamical problem. Indeed, this equation is rather complicated and hardly 
allows for a rigorous solution. Therefore, being interested in the qualitative understanding of 
statistics of the turbulent transport in the SOL, we develop a Monte Carlo discrete time simu- 
lation procedure which would help us to evaluate the asymptotic solutions of the given integral 
equation. 

General approach to the probability distributions of arrival times in such a discrete time 
model has been developed recently in [Hj. In general, its statistics can exhibit the multi- 
variant asymptotic behavior. Referring the reader to ^ for the details, in the present paper, 
we have just shown that the statistics of arrival times for the unstable fluctuations is either 
exponential or bounded by the exponentials (in particular, the latter would be true in the case 
of the randomly roaming wall) that is in a qualitative agreement with the data of numerical 
simulations and experiments |3]. 



2 Stochastic models of turbulent transport in the cross-filed system. 

The stochastic models of cross field turbulent transport used in the present paper refers to 
the effectively two-dimensional fiuid model of plasma based on the interchange instability in the 
SOL 13 ini recently discussed in '4]. In this model, one assumes the temperatures of ions and 
electrons to be constant, Tj <C Tg. Then the problem is reduced to that of two coupled fields, 
the fiuctuations of normalized particle density n{x,y,t) and that of vorticity field w{x,y,t), 
governed by the following equations 

Vtw = uoiyA±w - ^ — - dyu'' + /^, 

fc>i 

Vtn = uA^n + fn (2) 

written in the polar frame of reference with the normalized radial x = {r ~ a)/ p^ and y = 
aO/ ps poloidal coordinates. Time and space are normalized respectively to fij"^, the inverse 
ion cyclotron frequency, and to ps-, the hybrid Larmor radius. The covariant derivative is 
Vt = 9t + e^ ■ V X V, in which V = {dx-, dy), and A^ is the Laplace operator defined on 
the plane transversal to the axial magnetic field. The effective drive oc dy log {1 + n/ n) risen 
in the cross-field system due to the curvature of magnetic lines is represented by the series 
in dy n'' with the coefficients gk ~ ps / Rfif^ where R is the major radius of torus and n is 
the mean normalized particle density. The curvature coefficients gk averaged along the lines 
of magnetic field are considered to be constant and small parameters in the problem. The 
diffusion coefficients v and Uqv both are normalized to the Bohm's value Tg /eB and govern the 
damping in the small scales, herewith mq is the dimensionless Prandtl number where the knot 
distinguishes its value in the free theory from its effective value u in the renormalized theory 
(see Sec. ID). 

The Gaussian distributed random forces /„ and /^ in (j21) play the role of stochastic sources in 
the dynamical problem maintaining the system out of equilibrium and simultaneously modelling 
the Bohm's boundary conditions at the sheath which have not been explicitly included in Q, 
in contrast to the original models |l]-jS|- Herewith, the physically important effect of particle 
escape at the sheath is replaced with a quenched loss of particles in the SOL at the points 
for which fn{x,y, t) < 0. Simultaneously, the particles are supposed to arrive in the SOL 
in the areas where fn{x,y, t) > modelling the injection of particles from the divertor core 
along with the perturbations risen in the system due to the Langmuir probes 0-111] • For a 
simplicity, in the present paper, we assume that the processes of gain and loss of particles 
are balanced in average therefore (/„) = 0. The stochastic source of particles is used instead 



of the continuously acting radial Gaussian shaped source (localized at x = 0) studied in the 
numerical simulations j3] . Similarly, we impose the random helicity source /^ exerting onto the 
vorticity dynamics in (P)). 

Furthermore, the random sources /„ and f^ account for the internal noise risen due to the 
microscopic degrees of freedom eliminated from the phenomenological equations (j21) . From the 
technical point of view, the random forces help to construct a forthright statistical approach 
to the turbulent transport in the SOL. In particular, it allows for the quantum field theory 
formulation of the stochastic dynamical problem Q (based on the Martin-Siggia-Rose (MSR) 
formalism ^2|) that gives a key for the use of advanced analytical methods of modern critical 
phenomena theory |13j . 

The Gaussian statistics of random forces in Q is determined by their covariances, 

A.n(r - r', t - t') = (/„(r, t)/„(r', t')) , D^^{r -r',t- t') = (/^(r, t)/^(r', t')) , r = {x,y), 

describing the detailed microscopic properties of the stochastic dynamical system. In the 
present paper, we primarily discuss the large scale asymptotic behavior of the response func- 
tions ( 6n{r, t)/ Sfn{0, 0) ) and ( 6n{r, t)/ 6fw{0, 0) ) quantifying the reaction of system onto the 
external perturbation and corresponding to the r- distributions of particle density fluctuations 
expected at time t > in a response to the external disturbances of density and vorticity 
occurring at the origin at time t = 0. The high order response functions are related to the 
analogous multipoint distribution functions F„ (ri,ti, . . . r„,t„; r'i,t[ . . . r'„,t„) as 

5" [n(ri,ti) . . . n{rn,tn)] \ ■^^ , , / / \ 

^ ' ^ ' / / permutations 

with summation over all n! permutations of their arguments ri,ti . . . r„, t„. 

We consider a variety of microscopic models for the random forces /„ and fw in the stochastic 
problem Q. Under the statistically simplest "white noise" assumption, these random forces 
are uncorrelated in space and time, 

Dnn (r - r', t - t') = Tn 6{r' - r) 6{t - t'), D^^ (r - r', t - t') = T^ 6{r' - r) 6{t - t'), (3) 

in which r„ and F^ are the related Onsager coefficients. 

Recent studies reported on the statistics of transport events in the cross-field systems |lj,|Ilj 
pointed out the virtual importance of correlations existing between density and vorticity fiuc- 
tuations in the dynamical problem. In particular, this effect is referred to the formation of large 
density blobs of particles close to the divertor walls by attracting particles via the cross field 
fiow, the latter being the larger for strong blobs with strong potential gradients |3]. Indeed, in 



the physically realistic models of turbulent transport in the SOL, it seems natural to assume 
that the random perturbations enter into the system in a correlated way. To be specific, let us 
suppose that there exists a finite reciprocal correlation time Tc (| r' — r |) > between the ran- 
dom sources /^(r,rc) and /„(r',0) in the stochastic problem (j21). For a simplicity, we suppose 
that the relevant relaxation dynamics is given by the Langevin equation, 

'dt ~ n^T^' ^^^ 

in which /5 ~ ( /^ ) > 0. In the momentum representation, the non-local covariance operator 
T^^ can be specified by the pseudo-differential operator with the kernel 



r. 



-\k) = Xiyk^-^\ < 27 < 1, (5) 



which specifies the characteristic viscoelastic interactions between the "fast" modes of density 
and vorticity fiuctuations. The coupling constant A > naturally establishes the time scale 
separation between "fast" and "slow" modes. In the case of 27 <^ 1, the Langevin equation (JH) 
with the kernel (J3J reproduces the asymptotical dispersion relation typical for the Langmuir 
waves travelling in plasma, uj ~ A;^"''* as A; ^ with rj^, ~ 0.0804 (in three dimensional 
space) pS]. Alternatively, for the exponents 27 -^ 1, it corresponds to the ion-acoustic waves 
travelling in the coUisionless plasma with the velocity Xu ~ ^/Te / M where M is the ion mass. 
Intermediate values of 7 correspond to the various types of interactions between these two types 
of plasma waves described by the Zaharov's equations [T^ . 

The relaxation dynamics (j^-E)) establishes the relation between the covariances of random 
sources in (J2I), 

L'„„(r,t) = — — drdt / dp — — L)^^(r -r, t -t), r = |r|, (6) 

47r/5 J Jq Aup'- ^^ 

where Jq is the Bessel function of the first kind. In the present paper, we choose the covariance 
of random vorticity source, 

(/^(r,t)/^(r',t')) = [^ f -^D^^iu,k)exp[-tiu{t-t')+tkir-r')], k = \k\, 

in the form of white noise (jHI), in which the relevant Onsager coefficient F^ is found from the 
following physical reasons. Namely, the instantaneous spectral balance of particle fiux, 

W{k) = lJ^{Uk,u;)U-Ku;)), (7) 

derived from (jHI) should be independent from the reciprocal correlation time Tc{k) at any k 
that is true provided D^w {^, k) oc X k~^^. Furthermore, the Onsager coefficient T^ has to 



fit into the appropriate physical dimension which is assembled from the relevant dimensional 
parameters, Uqu and k. Collecting these factors, one obtains the ansatz D^^y^ oc Xu^u^ k^~'^~'^'^ , 
in which d = 2 is the dimension of space. The power law model for the covariance of random 
helicity force oc k^-'^-'^'y does not meet the white noise assumption since 6{r — r') ~ k^ and 
therefore calls for another control parameter 2e > 0. Eventually, in the present paper, we use 
the model 

D^^ik, uj) = T^k^~''-^'-^\ r^ oc Xuy, d = 2, (8) 

with the actual value of regularization parameter 2e = 4, for d = 2. 

Let us note that the ansatz (jH} is enough flexible to include the various particular models 
of particle pump into the SOL. For instance, the alternative to the white noise assumption 
spatially uniform particle pump for which the covariance 

in the large scales, can be represented by the ansatz (jH} with the actual value 2e = 3. 

In the rapid-change limit of the stochastic model, A ^ oo (i.e., Tc -^ 0), the covariance (jHl) 
turns into 

( /„(k, .;) /„(-k, ^) ) ^ ^ e-'''^^'\ (9) 

and recovers the white noise statistics Q along the line e = j, in d = 2. Alternatively, in the 
case of A ^ (that corresponds to Tc -^ oo), the time integration is effectively withdrawn 
from ©, so that the resulting configuration relevant to (jHl) appears to be static oc A;'^~'^-2e ^^-^^ 
uncorrelated in space (at d = 2) for 25 = 2. 

The power-law models for the covariances of random forces has been used in the statistical 
theory of turbulence J7] (see also the references therein). The models of random walks in 
random environment with long-range correlations based on the Langevin equation (0)) have 
been discussed in concern with the problem of anomalous scaling of a passive scalar advected 
by the synthetic compressible turbulent flow [l8j, then in ^\, for the purpose of establishing 
the time scale separation, in the models of self organized criticality 



3 Iterative solutions of the stochastic problem and their diagram 
representation 

The linearized homogeneous problem, for the fluctuations of density n and vorticity w van- 



ishing at t — i> oo, 

[dt - v/^^]%n = 5{v)5{t), [dt - u^v A^]%u,{r, t) + gidy'Knij, t) = 5{v)5{t), (10) 



is satisfied by the retarded Green's functions, 



A„(r,t) 






Aixut 



exp 



Avt 



A^(r,t) 



d{t) 

Anuut 



exp 



r 



Auut 



29ie{t) X 
u (u + 1) r^' 



fitting into the retarding conditions, A„(r,t) = A^(r,t) = 0, for t < 0, that express the 
casualty principle in the dynamical problem. Nonlinearities in (j21) can then be taken into 
account by the perturbation theory. 



n(r, t) 



J dv'dt' K(r-r', t-t') [/n(r',t') - v(r',t') x Vn{v\t') 



w{r,t) 



dr'dt' AJr~r\t-t') 



Ur',t' 



k>l 



'ir9k 



dyn\r\t')-w{Y\t') xVw{r',t' 



+gi I dv'dt' X^(r-r', t-t') f dr"dt" S:„(r'-r", t' -t") dy [/„(r", t") - v(r", t") x Vn(r", t")] . 

(11) 
The solutions (fTTj) allow for the diagram representation (see Fig. Q), where the external line 

(a tail) stands for the field n, the double external line denotes the field w, and the bold 

line represents the magnetic flux v. The triangles stay for the random force /„, and the filled 

triangles represent fw The retarded Green functions A„ are marked by the lines with an arrow 

which corresponds to the arguments (r', t') and (r", t") (the direction of arrows marks the time 

ordering of arguments in the lines). Similarly, the double lines with an arrow correspond to 

the retarded Green functions A^. Slashes correspond to the differential operator V. Circles 

surrounding vertices representing the antisymmetric interaction v x V, squares present the 

vertices proportional to the poloidal gradient dy. 

All correlation functions of fluctuating fields and functions expressing the system response 

for the external perturbations could be found by the multiplication of trees (jll|) displayed on 

Fig. ^ followed by the averaging over all possible configurations of random forces /„(r,t) and 

fw{r,t). In diagrams, this procedure corresponds to all possible contractions of lines ended 

with the identical triangles. Thereat, the diagrams having an odd number of external triangles 

(correspondent to the random forces) give zero contributions in average. As a result of these 

contractions, the following new elements (lines) appear in the diagrams of perturbation theory: 

A„„(r-r',t'-t) = / dridti / dr2dt2 A„(r - ri, ti - t)D„„(ri - r2, ti - ^2) A„(r2 - r', t2 - 1'), 



A^^(r-r',t'-t) = / dridti / rfrarfta A^(r -ri, ti -t)Z:'^^(ri -rg, ti -^2) A^(r2 -r', ta -i^') 

+gld^yA^n{r-r',t'-t), 
A^n{r - r', t' -t) = gi / dridti An(r - ri, ti - t)(9yA^^„(ri - r', ti - t'), 

^^„(r - r', t'-t)=gi dndh K(r - ri, ti - t)9^X(ri - r', h - t'), (12) 

which are the free propagators of particle density and vorticity fluctuations, the mixed corre- 
lator, and the retarded mixed Green's function. In diagrams, we present the free propagators 
(I12|) by the correspondent lines without an arrow, and the retarded mixed Green's functions 
by the composite directed lines (see Fig. |2I). 

The cross-field drift function v(r,t) is not involved into the linear homogeneous problem 
(fTUj) and, therefore, it does not appear in the free propagators (O, however, it is presented in 
the nonlinear part of dynamical equations and therefore appears in the diagrams of perturbation 
theory. Due to the simple relation w = V x v, the propagators containing the field v are 
the same as those with w: Aw, A^v, and A„v- The bold lines representing v in Fig. Qcan 
be replaced in the diagrams of perturbation theory with the double lines (which correspond to 
the field w) with the additional factor (in the momentum representation) —iSzms km/k"^-, where 
m,s = x,y and Ezxy is the antisymmetric pseudo-tensor, for each v. 

In this framework, the exact correlation functions of fields and the response functions can 
be found from the Dyson equations, 

i^S-K^-^. (|;i)-.K'-.. (13) 

where S„ and S^ are the infinite diagram series, in which the first diagrams are shown in Fig. |21 
The diagram technique introduced in the present section is suitable for the system preserving 
the continuous symmetry of (0), apart from the sheath. 

4 Functional integral formulation of the stochastic problem. Dimen- 
sional analysis, UV divergencies in diagrams and renormalization 

In the present section, we study the properties of diagram series resulting from the iterations 
of the stochastic dynamical equations with the consequent averaging with respect to all possible 
configurations of random forces. The diagrams for some correlation functions diverge in small 
scales. The use of conventional arguments borrowed from the quantum field renormalization 
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group jzd] helps to prove the consequent subtraction of the logarithmic divergent terms in all 
orders of perturbation theory out from the diagrams. 

The set of diagrams arisen in the perturbation theory by the iterations of (j21) is equivalent 
to the standard Feynman diagrams of some quantum field theory with the doubled set of 
stochastic fields: the fluctuations n and w, the flux function v, the auxiliary fields n', w' 
functionally conjugated to the Gaussian distributed random forces /„ and /u; in Q, and the 
Lagrange multiplier v' for the binding relation w = V x v. The coincidence of diagrams is 
a particular consequence of the general equivalence between the t-local stochastic dynamical 
problems (in which the interactions contain no time derivatives) and the relevant quantum 
field theories [25 with the action functional S found in accordance to the MSR formalism 
[T^ . Statistical averages with respect to all admissible configurations of random forces in a 
stochastic dynamical problem can be identified with the functional averages with the weight 
expiS. In particular, for the stochastic problem (j2)), the generating functional of the Green 
functions, Q{Aq>), with the arbitrary source fields A<j,(r, t) where $ = {n, n', w, w', v, v'} can be 
represented by the functional integral 



g{A^)= j V<^ exp S{^)+ I drdt A^{r,t)<l>{r,t) 



(14) 



in which 



5($) = - dr dr' dt dt' [n'{r, t)L'„„(r -r',t- t')n'{r', t') + u;'(r, t)D^^{Y - r', t - t')w'{r\ t')] 

(15) 
— / drdtn{r,t) [dtn{r,t) + v(r,t) x Vn(r,t) — h'A±n{r,t)] 



-ifgt 



dtw{r, t) + v(r, t) x Vw(r, t) + V ^ ^j ^^ dyn''{r, t) - UQuA^w{r, t) 

K 
k>l 



+ / drdtV{r,t)[w{r,t) - V x v(r,t)]. 

The source functions v4„/( r, t ) and A^i{r, t ) in (J14j) are interpreted as the not random external 
forces, so that the Green functions ( w{r, t) w'{r', t')) , { n(r, t) n'(r', t') ) , and ( n{r, t) w'{r', t') ) , 
for t' < t, coincide with the response functions {5w /5fw ), {6n/6fn ), and {5n jdfw ) re- 
spectively. All possible boundary conditions and the damping asymptotic conditions for the 
fluctuation fields n and w at t — i> — oo are included into the functional integration domain in 

The functional integral formulation ()14p of the stochastic dynamical problem Q allows 
for the use of various techniques developed in the quantum field theory to study the long- 
time large-scale asymptotic behavior of the quantum and stochastic systems. The integral 
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()14|) has the standard Feynman diagram representation which is equivalent to the iterative 
solution of © with an exception of the self-contracted lines Tr f A$$'j = Tr f A$/<j,j = 0, 
which are proportional to oc 9{t — t') in time-representation and therefore discontinuous at 
t = t'. Stipulating that 9{0) = 0, one can exclude all redundant graphs from the perturbation 
theory of ()14|) . Lines and vertices in the graphs of perturbation theory are defined by the 
conventional Feynman rules and correspond to the free propagators (equivalent to (fT^ ) readily 
calculated from the free (quadratic) part of functional (fT3j) and the nonlinear interactions of 
fields respectively. In the actual calculations, it is convenient to use the propagators (fT^ in 
their momentum-frequency representation, 

A„n'{k,Lj) = {-iuj + vk'^)"^ , A^^>{k,uj) = {-iu + uouk^y^, A^v'(^) = -£zmsks/k^, 

A.nn[k,UJ) = l\nn' J^nn\,kjOj) iXn'ri} Ay^yj^k^LO) = g^ky/\nn[k, OJ) + A^„^/iJ^^(A;, CUJ A^;,/^, 

Anwik,uj) = igiAnn'ik,uj)ky A^^,{k,uj), Anw'{k,uj) = igiAnn'{k,uj)ky A^^,i{k,uj), 

Ann'{k,uj) = A„„,(A;,u;), A^^^{k,uj) = A^^,{k,uj), A„^/(fc,u;) = A^,^{k,uj). 

Propagators including the field v coincide with those of the field w up to the multiplicative 
factor —iSzms km/k'^ for each field v. Propagators of auxiliary fields A^iij,i{k,uj) = 0. 

The action functional p5|) is invariant with respect to the generalized Galilean transforma- 
tions of fields in the poloidal direction, 

Vy{r,t) -^ Vy{r,t) - a{t), n{x,y,t) -> n{x,y - h{t),t) (16) 

with the parameter of transformations a( t ) (the integrable function decaying at t ^ — oo) and 
b(t) = J_ a(t') dt'. Furthermore, any quantity Q in ()15|) can be characterized with respect 
to the independent scale transformations, in space and time, by its momentum dimension dq 
and the frequency dimension dq. In the "logarithmic" theory (which is free of interactions that 
is analogous to the linearized problem (fTUj) ). these scale transformations are coupled due to the 
relation dt ~ V^ in the dynamical equations. 

This allows for the introduction of the relevant total canonical dimension dq = dq + 2dq 
and the analysis of UV divergencies arisen in the diagrams of perturbation theory based on 
the conventional dimension counting arguments [2HII2H1- In dynamical models, dq plays the 
same role as the ordinary (momentum) dimension in the critical static problems. Let us note 
that the poloidal gradient term oc dy in (fTUI) is responsible for the stationary contributions 
into the Green's function as t > 0, so that the above definition of dq remains unambiguous. 
Stipulating the natural normalization conventions, d^ = —d^ = 1, d'^ = d'^ = 0, d^ = d^ = 0, 
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d^ = —df = 1, one can find all relevant canonical dimensions from the simultaneous momentum 
and frequency scaling invariance of all terms in (jl5|) (see Tab. 1). 

Integrals correspondent to the diagrams of perturbation theory representing the 1-irreducible 
Green's functions G diverge at the large momenta (small scales) if the correspondent UV- 
divergence index 6g is a nonnegative integer in the logarithmic theory, 

6G = d + 2- {N^d^ + N^,d^>) > 0, (17) 

where d is the dimension of space, (i$ $/ are the total canonical dimensions of fields $ and $', 
and iV$^$/ are the numbers of relevant functional arguments in G. As a consequence of the 
casualty principle, in the dynamical models of MSR-type, all 1-irreducible Green's functions 
without the auxiliary fields $' vanish being proportional to 6{0) and therefore do not require 
counterterms ^7j. Furthermore, the dimensional parameters and external momenta occurring 
as the overall factors in graphs also reduce their degrees of divergence ()17|) . 

In spite of the following 1-irreducible Green functions could be superficially divergent at 
large momenta: {n'n\), {n'\), (n'vv), {n'n . . .n) (with an arbitrary number of fields n ), the 
only Galilean invariant Green's function admissible in the theory (jl5|) which actually diverges 
at large momenta is {n'n) . The inclusion of the relevant counterterm subtracting their superfi- 
cially divergent contribution is reproduced by the multiplicative renormalization of the Prandtl 
number 

Uq = uZu 

where Uq and u are the bare and renormalized values of Prandtl's number. In principle, the rel- 
evant renormalization constant Z^ can be calculated implicitly from the graphs of perturbation 
theory up to a finite part of the relevant counterterm. 

However, the standard approach of the critical phenomena theory is useless for the deter- 
mining of the large scale asymptotes in the problem in question, since the severe instability 
frustrates the critical behavior in the system preventing its approaching to the formal asymptote 
predicted by the conventional renormalization group method. As a consequence, the critical 
dimensions of fields and parameters which can be computed from the renormalization procedure 
would have just a formal meaning. 

5 Large scale instability of iterative solutions 

The iterative solutions for the stochastic problem (j2)) constructed in Sec. 01 would be asymp- 
totically stable in the large scales provided all small perturbations of both density and vorticity 
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damp out with time. In particular, the exact response functions found from the Dyson equations 
(fT^ should have poles located in the lower half-plane of the frequency space. 

The stability of free response function {6n/ 6fn) = {—iuj + z/Zc^)"^ which effectively cor- 
responds to the linearized problem (jlU|) is ensured in the large scales by the correct sign of the 
dissipation term uk^ > 0. In a "proper" perturbation theory, apart from a crossover region, 
the stability of exact response functions is also secured by the dissipation term which dominates 
over the dispersion relation in the large scales, 

u{k) = —iuk"^ + iTjnn{k,uj) (18) 

where the self-energy operator S„„(fc, cj) is the infinite series of all relevant 1— irreducible 
diagrams of perturbation theory. However, in the perturbation theory discussed in Sec. El the 
leading contribution into the self-energy operator is S„n ~ 0{ky) that could lift up the pole of 
the response function into the upper half-plane of the complex cj— plane rising the instability in 
the system as fc — > 0. Such an anomalously strong contribution comes from the diagrams which 
simultaneously include both the antisymmetric interaction vertex ~ z ez ■ Eijz Vikj together with 
the poloidal gradient ~ iQkky and the free propagator of particle density A„„(A;,co'). Such 
diagrams appear in all orders of perturbation theory for the response functions ( 5n/ 6fn ) 
and ( 6n/ 5fu, ) indicating that the instability could arise due to the random fluctuations of 
both density and vorticity. The first order 1-irreducible diagrams for them are displayed in 
Fig. El The small parameters Qn with n > 2 appear in the forthcoming orders of perturbation 
theory. In particular, the last diagrams in both series shown in Fig. El generate the "anomalous" 
contributions in the large scales. 

To be certain, let us consider the second diagram in the series for ( 5n/ 6fn ) which corre- 
sponds to the following analytical expression: 

VI r^ I. f ^P fdujf du' k^p, + kypy -pI-pI ^ 

^""li-W - 92ky J ^^2] ^ J ^ l^-p^2 Ann(P,^) A.(k- p, ^ - .; ). 

(19) 

Being interested in the 0{ ky )-contribution into ()19|) . one can neglect the fc— dependence in the 

integrand. The analytical properties of this contribution depends very much upon the certain 

assumption on the covariances of random forces since it changes the free propagator A„„. For 

instance, under the white noise assumption Q A„„ ~ (cj^ + z/^p^)~^, the integral in (fT^ 

diverges at the small scales (large momenta) for e < 1 and diverges at the large scales (small 

momenta) for e > 1. Introducing the relevant cut-off parameters, one obtains the anomalous 

contribution 

^nnli.ioop - S g^ ^_^ + e) (^ + 1) ' y^""! 
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in which u is the renormahzed value of Prandtl's number. In the preceding section, we have 
shown that the logarithmic divergencies risen in diagrams for the response function ( 6n/ 6fn ) 
in the small scales (large momenta) can be eliminated from the perturbation theory by the 
appropriate renormalization. The singularity in (j2Up arisen at the small momenta m ^ for 
e > 1 would compensate the smallness of g2, so that any density fluctuation with ky > 
appears to be unstable. 

Accounting for the finite reciprocal correlation time r^ > between the vorticity and density 
random sources in Q introduces the new dimensional parameter A > into the particle density 
propagator, A„„ ^ (u;^ + u'^p'^ (1 + \p~'^'^) ) . Then, the integral (fT^ can be computed by 
its analytic continuation for any momenta excepting for the isolated points, —l+e = 7 mod 1, 

The dispersion relation p8|) determines the region of asymptotic stability in the phase space of 
cross-field transport system. Namely, a density fluctuation arisen in the SOL with some random 
momenta {k^, ky) would be asymptotically stable with respect to the large scales fc — ;> if 

^' < ^w^ (22) 



kl + kl A(e,7) 

and be unstable otherwise. In the first order of perturbation theory, the amplitude factor A{e, 7) 
is given by (PT|) . For different values A{e, 7), the stability condition (j^^ determines the set 
of circles (see Fig. H)) osculating at the origin which bound the unstable segments of phase 
space. One can see that the density fluctuations with ky -^ (i.e. extended in the poloidal 
direction) are asymptotically stable for any \kx\ > 0. Density fluctuations characterized by 
\ky\ > would be asymptotically stable in a certain stochastic model provided sign (ky) = 
— sign {A{ e, 7 )) , for the given values of e and 7. The signature of the 1-loop order contribution 
()21|) into A{e,'y) is displayed on the diagram in Fig. El (black is for +1, white is for —1) at 
different values of 7 and d = 6 — 26 — 27, the space dimension related to the actual value of 
regularization parameter e under the statistical assumption (jH}. 

It is important to note that the stability condition ()22|) can be formulated as an upper bound 
for the order parameter ^ = ky / k"^: C, < C,c where C,c = ^i^, l)~^- For the uncorrelated 
statistics of random forces, in the stochastic dynamical problem, A{e,'j) (pH) diverges as /c — » 
and therefore ^c ■"* 0. 

6 Turbulence stabilization by the poloidal electric drift 

To promote the stochastic cross-field turbulent transport system Q from the instability to a 
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stable regime, it seems natural to frustrate the symmetry which breaks the Galilean invariancy 
in (j2)). This can be achieved by generating a constant uniform drift in the poloidal direction 
"^y ^ ^y ~ ^ (^"7 biasing the limiter surface, (f){x) -^ (f){x) + xV) that would eradicate those 
configurations with the trivial poloidal component of electric drift Vy = 0. In general, the 
relevant dispersion equation u {k,V) = could have many solutions Vk for A; ^ 1. Herewith, 
the turbulence stabilization is achieved for the drifts V from the intervals Vk-i < V < Vk for 
which Im(co') < 0. 

To be certain, let us consider the dispersion equation correspondent to the simplest response 
function {6n/ 6fn), in the 1-loop order. The leading contribution into the dispersion equation 
is given in the large scale region by the diagram p9p . Under the white noise assumption Q , the 
free propagator accounting for the uniform electric drift V is A„„ ~ ( cj^ + z/^/c^ + V^ky"^ ) 
Then, for \ky\ < 1, the dispersion relation reads as following, 

^ ' ' ^ u + 1 VStt (-l+£) ^ 8z/ (1 + n) r(2 + £) 

The latter relation shows that for any finite rra > there exists finite Vc < oo such that for 
any V > Vc one obtains Im(u;) < in ((221) , however, Vc ^ oo as m — > 0. 

In contrast to it, in the case of correlated statistics (|1|-E1 El, the dispersion relation is 
not singular for /c ^ excepting some particular values of e and 7, and the correspondent 
stabilizing electric drift, in the 1-loop order, equals to 

V = TT VU -^ — ^^ — - (24) 

(l + M)'+('"+')/'^ + M(2e + l)-27 (1 + 2m) sin 715/7 

In the range < 7 < 1/2, m > 0, this expression is singular at the points 5/7 G Z. 

7 Qualitative discrete time model of anomalous transport in the 
SOL 

Large scale instability developed in the cross-field model (j21) is related to the appearance and 
unbounded growth of fluctuations of particle density close to the wall. In accordance to the 
fluctuation-dissipation theorem, the fluctuations arisen in the stochastic dynamical system are 
related to its dissipative properties. In particular, the matrix of the exact response functions 
M( k, cij ) expressing the perturbations of fields n and w risen due to the random sources /„ and 
fu) determines the matrix of exact dynamical Green's functions G( k, cj ), 

M(k, cu) - M^(-k, -cj) = icjG(k, cj) (25) 
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where M^ is the transposed M. In the large scale limit k ^ 0, we take into account for the 
leading contributions into the self-energy operators in the elements of M, 

Rnn — [—iuj — u Alky +z/fc^) , Rn^ ~ (— ztu — uuA2ky +h'uk^) , (26) 

R^n — {-i^ +iyk'^) , Rww = {-i^ +iyuk'^) , 

in which Ai^2 are the amplitudes of the anomalous contributions competing with the dissipation 
oc 0{k'^) in the large scales. Fluctuations of particle density arisen in the model (j21) grow up 
unboundedly provided either Ai ky > k"^ or A2 ky > k"^ for the given values of e and 7. The 
correspondent advanced Green's functions appear to be analytic in the lower half-plane of the 
frequency space, 

T^-^nnKKt)- ^4 _ ^2 ^2 ' '^ « ■ G„^ ( k, ij _ - _ —, (27j 

being trivial for t > 0. 

For instance, let us consider the advanced Green's function G^„ which relates the density 
of particles 6n{ r, t ) in the fluctuations characterized with Ai ky > k"^ and arisen at the point 
r G f2 inside the divertor at time t with the particle density ^(r', t') of those achieved the 
divertor wall at some subsequent moment of time t' > t at the point r' G dVt: 

Hr',t')\,^=[ dt [ drdnir,t)G:^ir'-r) {2i 

where 



t<t' 



Gl 



2uAi 



' ■ My'-y) J (Ai\r'-r\\ Ai{y'-y) f Ai\r' - r\ 
sm Jo ^^+ cos Ho< 



(29) 
in which |r — r'| = \/{x — x'Y + {y — y')"^, Jo and Hq are the Bessel and Struve functions 
respectively. The integral in the r.h.s. of ()28|) is finite for any Q provided | ^i | > 1, but for the 
compact fi as I Ai I < 1. To be specific, let us consider the circle Cji of radius R as the relevant 
domain boundary and suppose for a simplicity that the density of particles incorporated into the 
growing fluctuations inside the domain is independent of time and maintained at the stationary 
rate ^o( r)- Then the r— integral in the r.h.s of (j^Hj) can be calculated at least numerically and 
gives the growth rate B{R) for those density fluctuations, 

Sn{R,T) = T -BiR), (30) 

where r is the travelling time of the density blob to achieve the divertor wall that can be 
effectively considered as a random quantity. It is the distribution of such wandering times 

17 



that determines the anomalous transport statistics described by the flux pdf in our simplified 
model. The discrete time model we discuss below is similar to the toy model of systems close to 
a threshold of instability studied in [SJ recently. Despite its obvious simplicity (the convection 
of a high density blob of particles by the turbulent fiow of the cross field system is substituted 
by the discrete time one- dimensional (in the radial direction) random walks characterized with 
some given distribution function), its exhibits a surprising qualitative similarity to the actual 
fiux driven anomalous transport events reported in j3]. 

We specify the random radial coordinate of a growing fiuctuation by the real number x G 
[0, 1]. Another real number R G [0, 1] is for the coordinate of wall. The fiuctuation is supposed 
to be convected by the turbulent fiow and grown as long as x < R and is destroyed otherwise 
{x > R). We consider X as a random variable distributed with respect to some given probability 
distribution function P{x < u} = F{u). It is natural to consider the coordinate of wall R as 
a fixed number, nevertheless, we discuss here a more general case when R is also considered 
as a random variable distributed over the unit interval with respect to another probability 
distribution function (pdf) P{-R < u} = Q{u). In general, F and Q are two arbitrary left- 
continuous increasing functions satisfying the normalization conditions F{0) = Q{0) = 0, 
F(oo) = g(oo) = 1. 

Given a fixed real number rj G [0,1], we define a discrete time random process in the 
following way. At time t = 0, the variable x is chosen with respect to pdf F, and R is chosen 
with respect to pdf Q. li x < R, the process continues and goes to time t = 1. Otherwise, 
provided x > R, the process is eliminated. At time t > 1, the following events happen: 

i) with probability r], the random variable x is chosen with pdf F, but the threshold R keeps 
the value it had at time t — 1. Otherwise, 

a) with probability 1 — r], the random variable x is chosen with pdf F, and R is chosen 
with pdf Q. 

If X > R, the process ends; ii x < R, the process continues and goes to time t + 1. 

Eventually, at some time step r, when the coordinate of the blob, x, drops "beyond" R, the 
process stops, and the integer value r resulted from such a random process limits the duration 
of convectional phase. The new blob then arises within the domain, and the simulation process 
starts again. 

While studying the above model, we are interested in the distribution of durations of con- 
vection phases P-qiT] F,Q) (denoted as P{t) in the what following) provided the probability 
distributions F and Q are known, and the control parameter rj is fixed. The motionless wall 
corresponds to ?7 = 0. Alternatively, the position of wall is randomly changed at each time step 
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as r] = 1. 

The proposed model resembles to the coherent-noise models |2S]-|2Z1 discussed in connection 
with a standard sandpile model [201 in self-organized criticality, where the statistics of avalanche 
sizes and durations take power law forms. 

We introduce the generating function of P( r ) such that 



P(.) = 5^.-P(r), P{r) ^ ^^^^' 

T=0 



r! ds' 



(31) 



s=0 



(32) 



and define the following auxiliary functions 

/■oo 

K{n) = / F{u)''dQ{u), 5K{n) = K{n) - K{n + l), 
Jo 

p{l) = rf K{1 + 1) , for / > 1 , p(0) = , 

q(l) = (1 - i]f Kiiy--^ , for / > 1 , g(0) = , 

r(/) = r^' [r/ 5/r (/ + 1) + (1 - r^) fs:(/ + 1) 5R:(0)] , for / > 1, r(0) = , 

p = ri bK{\) + (1 - 7]) K{\) 5/^(0) . 
Then we find 

P{s) = 6K{0) +ps + ,f-— [f{s) + pp{s)qis) + p K{l)q{s) + K(l) q{s)f{s)] , (33) 

l-p{s)q{s) 

where p{s), q{s), f(s) are the generating functions corresponding to p{l), q{l), r{l), respectively. 
In the marginal cases t] = and rj = 1, the probability P{t) can be readily calculated, 

P^=o(r) = K{iy 6K{0), P^=i{t) = 6K{t). (34) 

The above equation shows that in the case of r/ = 0, for any choice of the pdf F and Q, the 
probability P{t) decays exponentially. In the opposite case r] = 1, many different types of 
behavior are possible, depending upon the particular choice of F and Q. 

To estimate the upper and lower bounds for P(r) for any rj, one can use the fact that 

if(l)" < K{n) < K{1) and < 6K{n) < K{1) , n e N. 

Then the upper bound for K{n) is trivial, since < F{u) < 1 for any u G [0, 1]. The upper 
bound for K{n) exists if the interval of the random variable u is bounded and therefore can be 
mapped onto [0, 1] (as a consequence of Jensen's inequality, and of the fact that the function 
u :-^ m" is convex on the interval ]0, 1[ for any integer n). The calculation given in ^ allows 
for the following estimation for the upper bound, 

P„(r) < r]" ^K{t) + (1 - 7])K{1) 6K{0) [r/ + (1 - r/) K{1)Y-' (35) 
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and, for the lower bound, 

Pvir) > rf 6K{t) + (1 - r/) K{iy 6K{0) (36) 

= r/-P,=i(r) + (l-77)P,=o(r). 

We thus see that, for any < 77 < 1, the decay of distribution P(r) is bounded by exponentials. 
Furthermore, the bounds (j35p and (|36|) turns into exact equalities, in the marginal cases t] = 
and 1] = 1. 

The simpler and explicit expressions can be given for P{t) provided the densities are uniform 
dF{u) = dQ{u) = du for all u G [0, 1]. Then the equations (jH^ give, 

P,.„(r) = 2-(-') P.,.,(.) = ^_J^_^. (37) 

For the intermediate values of rj, the upper and lower bounds are 

"^^ ' ^l-r])2-(- + i) < P(r)<i f^y. (38) 



(r + l)(r + 2) ^ " - ^ ^ - 2 V 2 

The above results are displayed in Fig. El 

8 Conclusion 

In the present paper, we have considered the stochastic model of turbulent transport in 
the SOL of tokamak. This problem allows for the concurrent symmetries, and the system 
exhibits a sever instability with respect to both perturbations either particle density or vorticity. 
Instability would reveal itself in the appearance of high density blobs of particles hitting into 
the reactor wall. 

The accounting for the dissipation processes introduces the order parameter ^ = \ky\ / {kl + A;^) 
and its critical value ^c such that the particle density fluctuation 5n{^) grows unboundedly with 
time as ^ > ^c and damps out otherwise. In the present paper, we compute the value of ^c, 
in the first order of perturbation theory developed with respect to the small parameter ps / Rn 
where ps is the Larmor radius, R is the major radius of torus, and n is the mean normalized 
density of particles. 

Our results demonstrate convincingly that the possible correlations between density and 
vorticity fluctuations would drastically change the value ^c modifying the stability of model. 
Characterizing the possible reciprocal correlations between the density and vorticity fluctua- 
tions by the specific correlation time Tc-, we demonstrate that any fluctuation of particle density 
grows up with time in the large scale limit {k — > 0) as r^ ^ 00 (the density and vorticity 
fluctuations are uncorrelated) and therefore ^c = 0. Alternatively, S,c> ^ provided t^ < 00. 
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The reciprocal correlations between the fluctuations in the divertor is of vital importance 
for a possibility to stabilize the turbulent cross field system, in the large scales, by biasing the 
limiter surface discussed in the literature before 0|. Namely, if .^c > 0, there would be a number 
of intervals [Vfc_i, V^] for the uniform electric poloidal drifts V such that all fluctuations arisen 
in the system are damped out fast. In particular, in the first order of perturbation theory, there 
exists one threshold value Vc such that the instability in the system is bent down a.s V > Vc- 
However, V^ ^ oo as ^c ^ 0. 

To get an insight into the statistics of growing fluctuations of particle density that appear 
as high-density blobs of particles close to the reactor wall, we note that their growth rates are 
determined by the advanced Green's functions analytical in the lower half plain of the frequency 
space. We replace the rather complicated dynamical process of creation and convection of 
growing density fluctuations by the turbulent flow with the problem of discrete time random 
walks concluding at a boundary. Such a substitution can be naturally interpreted as a Monte 
Carlo simulation procedure for the particle flux. Herewith, the wandering time spectra which 
determine the pdf of the particle flux in such a toy model are either exponential or bounded by 
the exponential from above. This observation is in a qualitative agreement with the numerical 
data reported in [T. 
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Figure 1: The diagrammatic representation for Eq. Hll() . 
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Figure 2: First diagrams for the simplest response functions {6n/5fn) and {6w/dfu 
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Figure 3: The simplest 1-irreducible diagrams eontributing into the self-energy eorreetions E„„ and S„tu for the 
response functions {dn/5fn) and {Sn/Sf^). The field indices denote the type of propagators and simultaneously 
the type of vertexes. The slashes mark the positions of derivatives V, and the skewed triangles denotes the 
inverse operator curl" . 
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Figure 4: . The admissibility condition (|22(l defines the set of elliptic curves which bound the unstable segments 
in the phase space. Density fluctuations with ky -^ 0, (i.e. extended in the poloidal direction) appear to be 
stable in the large scales for any | fc^; | > 0. Those fluctuations characterized by | fcy | > would also be stable 
in the large scales provided sign (ky) — — sign ( A{e, 7) ) for the given values e and 7. 
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Figure 5: . The signature of 1-loop order contribution H21|) into A{e,j) (black is for +1, white is for —1) at 
different values of 7 and d = 6 — 2e — 2j, the space dimension related to the actual value of regularization 
parameter e under the statistical assumption (|5]l. 
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Figure 6: . The distributions of wandering times near the wall in the discrete time model of Sec.|3 in the case 
of the uniform densities dF{u) = dG{u) = du for all u G [0, oo) at different values of control parameter ry. 
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